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Abstract 

We present the results for a model calculation of resonant two-magnon Ra- 
man scattering in a spin density wave (SDW) antiferromagnet. The resonant 
enhancement of the two-magnon intensity is obtained from a microscopic anal- 
ysis of the photon-magnon coupling vertex. By combining magnon-magnon 
interactions with 'triple resonance' phenomena in the vertex function the re- 
sulting intensity line shape is found to closely resemble the measured two- 
magnon Raman signal in antiferromagnetic cuprates. Both, resonant and 
non-resonant Raman scattering are discussed for the SDW antiferromagnet 
and a comparison is made to the conventional Loudon-Fleury theory of two- 
magnon light scattering. 

PACS numbers: 75.30.Ds, 75.30.Fv, 76.50.+g 

I. INTRODUCTION 

Soon after the discovery of high-Tc superconductivity |T| Raman scattering experi- 
ments were performed on the antiferromagnetic (AF) parent compounds La2Cu04 
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and YBa2Cu^0Q The analysis of the two-magnon Raman intensity has proven since 

to be a valuable tool for probing the collective magnetic excitations in these layered mate- 
rials. Common to all antiferromagnetic cuprates is a well defined two-magnon peak in the 
Raman intensity in Big and a weaker but still significant signal in Aig scattering geometry 
at a transferred photon frequency near 3000 cm~^. The frequency of the two-magnon peak 
has allowed an estimate for the unrenormalized AF exchange coupling between the copper 
spins in the Cu02 planes of about 136meV^ in La2CuOi, consistent with results of neutron 
scattering experiments for the spin wave velocity [0 and for the zone boundary magnon 
energy Pfl. Most of the magnetic properties of AF cuprates are well described by model- 
ling the undoped Cu02 layers by a spin ^ Heisenberg model on a square lattice. Yet, some 
anomalous features of the two magnon intensity profile have remained a puzzle: both, the 
asymmetric and broad lineshape in Big and the appearance of a two-magnon signal also in 
Aig and B2g geometry cannot be obtained within the traditional Loudon-Fleury theory for 



two-magnon Raman scattering |TO|-p!^]. 



Furthermore, two-magnon light scattering in AF cuprates is a resonant phenomenon and 
the scattering intensity as well as the line shape depend on the incoming photon frequency. 
The Loudon-Fleury theory is in principle a theory for non-resonant Raman scattering, as- 
suming a phenomenological coupling of the incoming and scattered photons to the localized 
spins of the antiferromagnet as described by the coupling Hamiltonian 

Hl_f = (^inc ■ Uij){Esc ■ Uij){Si ■ Sj) (l) 

<i,j> 

where Einc and Esc are the electric field vectors for the incoming and scattered photons, and 
Uij is a unit vector connecting spin sites i and j |T0|-p!2|. However, light scattering experi- 
ments on AF cuprate compounds so far have been performed with laser photon frequencies 
comparable to the charge transfer energy gap of these insulating materials |T3| , p^ . Therefore, 
photon induced transitions across the insulating energy gap are the natural origin for the 
resonant features of the two-magnon signal. A successful theory for resonant two-magnon 
light scattering must for this reason retain the charge degrees of freedom of the electrons. 



Along similar lines as in the recent work of Chubukov and Frenkel |T5[ we perform a model 
calculation for a spin density wave AF which allows to explore the resonant enhancement of 
the two-magnon Raman intensity. We calculate the scattering intensity using a microscopic 
description for the photon-electron coupling and for the creation of a magnon pair. Final 
state magnon interactions are included within a diagrammatic formulation based on the half- 
filled single band Hubbard model. We show that in this framework the frequency dependence 
of the photon-magnon vertices gives rise to an enhancement of the high energy part of the 
two-magnon spectrum. Several experimental features are explained as a consequence of the 
interplay between the two-magnon peak and resonance phenomena of the photon-magnon 
vertex function. 

The paper is organized as follows: In chapter II we start with a brief review of the 
spin density wave formalism for the Hubbard model at half-filling. In chapter HI and IV a 
detailed description of our diagrammatic approach is given and the basic coupling vertices are 
calculated. The extension to including final state magnon-magnon interactions is presented 
in chapter V. In chapter VI we evaluate results for the non-resonant case and compare them 
to the conventional Loudon-Fleury theory. In chapter VII we explore the experimentally 
relevant resonant case and calculate the scattering intensity with resonant electron-magnon 
vertex functions. A discussion of our results and the comparison to the experimental data 
are presented in chapter VIII. 

II. THE SPIN DENSITY WAVE STATE 

We start from the single band Hubbard model on a square lattice which is assumed to 
describe the low energy physics of the Cu02 layers. The Hamiltonian for the Hubbard model 
in standard notation is given by 

H = y e(k)4 cr ctct . . c. . (2) 

k,cr k,l,q 
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Here c^^^ destroys (creates) an electron with momentum k and spin a =],[. is the number 



of lattice sites and U is the on-site Coulomb repulsion. The tight binding dispersion of the 
square lattice with nearest neighbour hopping only is given by e{k) = — 4t7^, where t is the 
hopping amplitude and 7g = [cos^ak^) + cos{aky)]/2. Throughout the rest of the paper t 
and the lattice constant a are set to unity. 

At half-filling the nesting property e{k) = —e{k + Q) with Q = (vr, vr) leads to an 
instability of the Fermi sea of non-interacting electrons towards a commensurate spin density 
wave ground state with wave vector Q. Following the standard procedure as originally 
outlined by Schrieffer, Wen, and Zhang ||16| we introduce the staggered magnetization M in 
the SDW ground state {ip) by: 



M 



(3) 



to linearize the interaction part of the Hamiltonian 
tonian is diagonalized by the linear transformation 



The resulting Hartree-Fock Hamil- 



( 



7 



ka 



U(fc.(T) 



U(fc,(T) 



I 



which leads to 



(4) 



(5) 



Here 7^^^^ 7^^^'' destroy (create) quasi particles in the SDW conduction and valence band, 
respectively. The primed summation is restricted to momenta in the magnetic Brillouin 
zone (MBZ), i.e. to the momenta of the occupied Fermi sea for the non interacting system 
(e(A;) < 0), and the quasi particle energy dispersion is given by 



E{k) = ^Je^{k) + A^ 



(6) 



The SDW energy gap A = UM between the valence and conduction band is determined 
from the gap equation 



and the transformation amplitudes in (^) are given by = y|[1 + e{k) / E{k)] and = 

Due to magnetic umklapp scattering from the periodic SDW potential the single particle 
propagator is no longer momentum diagonal and it is conveniently expressed cls cl 2 X 2 
matrix with respect to the momenta q G MBZ and q + Q. In this notation the Hartree Fock 
c-particle propagator in the SDW state is written as 



' uj + e{k) a A > 



V 



A 



1 



(8) 



LU — e{k) ^ 

Alternatively we will in subsequent chapters also use the diagonal propagator matrix for the 
SDW conduction and valence band quasi particles. The transformation between the c and 
7 representation for the propagator matrices reads 



Gl{k] u) = V^{k, a)G^{k, u;)V{k, a) 



(9) 



with 



G^{k;uj) 



V 



g;j 



uj - E{k) + i6 
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(10) 



uj + E{k)-i6 I 

The collective spin wave excitations in the SDW state are determined by the poles of the 
frequency Fourier transform of the transverse susceptibilities 



= ^(V^|T^|(t)5Zf (0)|^) 



(11) 



where here cr = ±. The local spin raising and lowering operators = 5*^ + oiS^^ are 
represented in terms of fermion operators by 



a/3 



(12) 



where denotes the Pauli matrices with ^ = x,y. In the following we denote by x(^ 2x2 
matrices as in (1^). 



For the calculation of the dynamic transverse susceptibility matrix we account for the 
residual interactions between the quasi particles beyond the mean field approximation by 
summing the standard (RPA) ladder diagram series |T^. This leads to the matrix equation: 



cr, — cr / 
I 



with 



Xo ( 



k 



m 



k,k+q 



^HS+q'^kS+q 



Y ^^k,k+q^k,k+q 



k,k+q 



X 



± 



,uj + E{k) + E{k + q)-i5 a; - E{k) - E{k + q) + i6\ 



and the coherence factors 



(13) 



(14) 



h,k+q ~ '^k^k+q + '^k^k+q ' ^k,k+q 



U 



m+q + n^ 



k+q 



(15) 



In Eq. (|T^) the upper (lower) sign refers to the off- diagonal (diagonal) matrix elements, re- 
spectively. The spin wave dispersion Ugyj^q) follows from the condition det[l — f/x+~(g; u)] = 
0. In the strong coupling limit U ^ t the RPA susceptibility matrix takes a very transpar- 
ent form displaying explicitly the propagating spin wave excitations. An expansion of Xrpa 



up to the second order int/U and u /U leads to [|T7 



x:r{q-,uj) 



-2 J (1 — 7g-) aoj 

auo -2J(l+7g-) 



(16) 



In this strong coupling limit the spin wave dispersion uJswio) = 2Jy 1 — (75-)^ is identical 
to the linear spin wave (LSW) theory result of the spin | Heisenberg antiferromagnet with 
exchange coupling J = At'^/U ||16|-|T9|. 

In order to describe the coupling of the photons to the electrons we consider the Hubbard 
Hamiltonian in the presence of a weak transverse electromagnetic field. The Coulomb inter- 
action of the Hamiltonian remains unchanged but the vector potential A{f, t) of the photon 
field introduces a phase factor exp(i // A(r, t) ■ dr) into the kinetic energy. Expanding the 
kinetic energy part of the Hamiltonian up to second order in A yields [pO[] : 
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k,(7 <? 

1 

+ ^f^N ^ M-QiKpiqi + q2)Api-q2). (17) 

Here, we have introduced the current density operator j{q) with components 

. 1 de(k) : 

and the effective density operator 

In second quantization the components Aa of the vector potential are expressed in terms of 
photon annihilation {a^) and creation (ai) operators by 



^""^ ' (20) 



e is the polarisation unit vector of the photon and ujq- = cq is the photon frequency. Since 
the photon wavelength for visible light frequencies is large compared to the crystal lattice 
spacing and to all microscopic length scales of the electronic system we can safely neglect the 
photon momenta and henceforth use the q = limit for the photon field. The polarisation 
vectors of the commonly used experimental scattering geometries with linearly polarized 
light are collected in Table 1, where Cj and e/ denote the polarisation unit vectors of the 
incoming and scattered photon, respectively, with regard to the bonds of the CUO2 lattice. 

The Hamiltonian ( p!7D describes how the photon field couples to the current density 
operator and determines the bare coupling vertices for the common scattering geometries. 
The second order photon-electron coupling to the effective density will be neglected in the 
following, because it contributes to the Raman scattering intensity only at high frequencies 
due to direct interband transitions and does not lead to any resonance phenomena. 
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Symmetry 




J 


^ '■Lg 






Big 






B2g 


(1,0) 


(0,1) 



Table 1. Polarisation vectors for the incoming (cj) and 
scattered photon (e*/) in the common scattering geometries. 



III. TWO-MAGNON RAMAN SCATTERING 

As an alternative to the Golden Rule analysis of the Raman scattering intensity, we 
follow here a convenient diagrammatic formulation This method has been applied 



previously for calculating the Raman intensity from two-phonon scattering |]22| or for two 
spin fluctuation scattering in the paramagnetic phase of the Hubbard model [^. Here we 
extend this technique to a 2 x 2 matrix formulation for the SDW state. The general diagram 
for the Raman amplitude in the absence of magnon-magnon interactions is shown in Fig. 
|l]. The photons couple to the electronic system through a vertex function V , creating 
two magnons with momenta q and —q. The left side of the diagram corresponds to the 
physical process; the scattering intensity is then obtained from the symmetrically completed 
diagram by taking diagram cuts (see below) - a technique which yields results equivalent 
to the Golden Rule analysis. The virtue of this diagrammatic technique is that it allows 
to select those scattering processes which dominantly contribute to the two-magnon Raman 
signal p2[. Here and in the following we use the finite temperature Matsubara formalism 



and after analytic continuation to the real frequency axis we evaluate our results in the zero 
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temperature limit. 



-q 

FIG. 1. General diagram for the Raman amplitude of two-magnon scattering. The dashed 
lines represent the incoming and outgoing photons with frequencies iLO^,iu)y and the wiggly lines 
the magnon propagators. The vertices V and V contain the microscopic coupling of the photons 
to the magnon excitations and iuot is the transferred photon frequency. 

The final state of interest contains two magnons with opposite momenta. Therefore the 
two-magnon scattering intensity can be deduced by taking a cut of the diagram for the 
Raman amplitude across the two magnon lines only. After analytic continuation of the 
photon Matsubara frequencies the cut translates into taking the discontinuity of Fq across 
the real frequency axis with respect to the transferred photon frequency iujt-, 

J(Acj) = [Fo(ic<;t ^ At<; + i(5) - Fo(2cjj ^ Atu , (21) 

and /(Acu) is directly proportional to the scattering intensity. Explicitly, given the magnon 
propagators as represented by the transverse dynamic spin susceptibilities, and the vertex 
functions V and V\ which will be calculated in chapter IV, the Raman amplitude Fq has 
the following form: 

X'^'^i-Q, -q', i^t - iuJp)V^{q\ iuJt, iuJu, i^p) (22) 

where (3 is the inverse temperature. In order to single out only the two-magnon contribution 
to Fq, the residues of the vertex function must be disregarded when carrying out the internal 
frequency sum. Similarly, the cut prescription Eq. ( ^I]) with respect to iut is restricted to 



the arguments of the susceptibilities. (For details of this technique see Ref. |2^.) The 
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omitted high energy contributions correspond to final states with particle-hole excitations 
only. 

The frequency sum in ( p2D is conveniently performed by introducing the spectral repre- 
sentation for the transverse dynamic spin susceptibility 

X ' g;«^p = — / : du . (23) 

7r J lujp — uj 

The analytic continuation of the Matsubara frequencies in the vertex function is performed 



according to the rules i2^ . |23 



iuj^ - 


uJi + id 


iujy - 


—oji + i5 


iujt - 


UJi-UJf = 



J . (24) 

uji and ujf are the frequencies of the incoming and outgoing photon, respectively, and Acj = 
uji — ujf is the photon frequency shift. Taking the cut according to Eq. (|21| ) then yields the 
desired two-magnon contribution to the scattering intensity: 

I{Auj) = J duV„{q, Au, Ui + -u)\\axRPAi.q^ ^ X 

Imx^P^(-g, -q\ Au + ij)Vl{q\ Au, -uJi + iS, -uj) {n{uj + Au) - n{uj)} . (25) 



Here, the Bose distribution function n{u)) is evaluated in the T —>■ limit. Note that 
and are in fact the same vertex functions but differ in the argument for the incoming 
photon frequency due to the analytic continuation rules Eq. (|^). As a consequence the 
primed vertex function V^{q) is equal to the complex conjugate of Va{q). In order to outhne 
the derivation of the scattering intensity (^) no final state magnon-magnon interactions 
have been included so far. The inclusion of magnon-magnon interactions will be discussed 
in chapter V. 
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IV. THE VERTEX FUNCTION 



Now we focus on the photon-magnon vertex function which depends on the incoming 
hght frequency and is thus responsible for the resonant behavior of the scattering intensity. 
The electron-photon coupling Hamiltonian (|1^ determines the bare coupling vertices of the 
photon's vector potential to the current density. Given the bare photon-electron vertices the 
simplest and in the strong coupling limit most relevant contributions to the vertex function 
Vo-, which mediates the indirect coupling of the photon to the spin wave excitations, are 
shown by the diagrams in Fig. 0. 




FIG. 2. a.) Effective vertex function for the photon-magnon coupling. The solid, dashed, and 
wiggly lines denote the SDW c-fermion, the photon, and the magnon propagator, respectively. The 
spin labeling of the fermion lines has been omitted, it is explicitly indicated in b.). The solid circle 
represents the bare photon-electron coupling j- A, the filled triangle represents the electron-magnon 
coupling where the magnons are contained in the particle- hole ladder series as shown in b.). 

The algebraic expressions corresponding to the individual diagrams that contribute to the 
photon-magnon vertex function as shown in Fig. ^ are explicitly listed in the appendix. 
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Since the total vertex function is independent of a we henceforth drop the spin index 
cr. In addition to the diagrams shown in Fig. |^a the diagrams with reversed direction 
of the fermion hnes have to be included as well. With Eq. ( pSj) and the vertex function 
contributions as listed in the appendix the derivation of the Raman intensity in the absence 
of magnon-magnon interactions is complete. 



V. MAGNON-MAGNON INTERACTION 

The remaining step is to include the effects of final state magnon-magnon interactions. 
From the results obtained in the framework of the Heisenberg model it is well known that it 
is crucial to include final state interactions in order to avoid a diverging Raman intensity at 
twice the maximum magnon frequency Auo = 2u!^^ = 4 J resulting from the square root di- 
vergence of the magnon density of states at the MBZ boundary [|T0]-[12|. We therefore extend 
the Raman amplitude Fq to allow for repeated magnon-magnon scattering. This is achieved 
by replacing a bare vertex V in the Raman amplitude by a renormalized magnon-photon 
vertex function F which contains an infinite series of magnon pair scattering processes. Di- 
agrammatically this is represented in Fig. |]a and requires the solution of a Bethe-Salpeter 
like equation for F (see Fig. [po). The diagrams for the Raman amplitude F and the vertex 
function F translate into the following equations: 

Xrpa(-'?i> -Qi, - i^p)T{qi, iuj^, iuJt, i^p) , (26) 



r{q[, iUf,, iuJt, iuJp) = V{q{, iu^, iut, iujp) -j^Y^^Yl ^2, ^^p, iuJt, iuJi) x 

^ ii^l 92,'? 2 

Xrpa{Q2, qL i^i)XRPA{-Q2, -qL i^t - iuJi)T{q2, iuj^, iuJt, iuJi) . (27) 
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a.) 



F = 




b.) 






+ 



FIG. 3. a.) General diagram for the Raman amplitude with final state magnon-magnon in- 
teraction, b.) Diagram for the renormalized magnon-photon vertex function T which contains 
only magnon-number conserving scattering processes. The solid squares in a.) and b.) repre- 
sent the bare photon-magnon coupling V (see Fig. |^) and V'^ in b.) represents the irreducible 
magnon-magnon interaction. 

For the irreducible magnon-magnon interaction we include only magnon number 
conserving scattering processes which - at least for the Loudon-Fleury theory of Raman 
scattering from the Heisenberg antiferromagnet - have been shown to be the most important 



2J]. For our present SDW state based calculation the magnon-magnon interaction vertex 



has to be expressed in terms of the residual Hubbard interaction between the fermionic 
quasi particles. The simplest diagrams which serve for this purpose are shown in Fig. ^ and 
involve an internal loop with 4 fermion Green's functions. Fortunately, for our subsequent 
strong coupling evaluation of the two-magnon Raman intensity in Big geometry we have 
indeed found that these two diagrams give the dominant contribution to the magnon-magnon 
interaction. As we will show below, evaluating in the static zero frequency limit the result 
for the Big intensity can still be obtained analytically for U ^ t in an intermediate frequency 
range despite the mathematical complexity of the Raman amplitude F in Eq. (^6|). 
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^1 



V(^J,) = 



-q 



1 '^2 M ^2 

FIG. 4. Vertices for magnon-magnon scattering. As before (see Fig. ^a) the filled triangle 
corresponds to the electron-magnon coupling. 



VI. NON-RESONANT SCATTERING 



First we analyse the results in the non-resonant case, when the incoming photon fre- 
quency uji is much smaller than the SDW energy gap 2 A. In the strong coupling limit all 
expressions take a transparent form and become analytically tractable. In Big symmetry 
the rhomb-shaped vertex diagram A in Fig. |^ gives the only contribution to the total vertex 
function; the diagrams -B,C, and D vanish identically to leading order in t/U. Explicitly 
we find 



(28) 



with 



[cos(g^) - cos( 



A2 

and /(.;,) oc _ 



(29) 



This simple structure of the photon-magnon vertex function leads immediately to a simple 
expression for the scattering intensity (^) in Big geometry without final state magnon- 
magnon interactions: 



[UJt 



4A 



2^2 



^ 1 



7o- 



-6{Auj - 2uJswiq)) 



(30) 



As expected the neglect of magnon-magnon interactions results in a logarithmic divergence 
of the Big intensity at twice the maximum spin wave frequency 4J. 

In Aig symmetry all 4 diagrams in Fig. |^ for the vertex function give non-vanishing 
contributions, but due to a perfect cancellation of the different terms the total vertex function 
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and therefore the scattering intensity vanish to leading order in t/U . This cancellation has 
been previously noted by Chubukov and Frenkel |[T^. Also the intensity in symmetry 



is zero in the strong coupling limit. For weak and intermediate values of f//t, however, the 
two-magnon intensity is finite in all three symmetry channels. 

The logarithmically diverging scattering intensity in Big ( PD| ) and the vanishing inten- 
sities in Aig and B2g symmetries for U ^ t are results one also obtains using the effective 
Loudon-Fleury spin-photon coupling Hamiltonian in the framework of the spin | Heisenberg 
model. In the non-resonant region, where (2A— cuj) ^ J, and in the strong coupling limit the 
assumption of localized spins is valid and the Loudon-Fleury approach is expected to yield 



the correct description |1T5| , P0| . Since we have demonstrated the equivalence between the two 
approaches at least in the absence of magnon-magnon interactions, the correct selection of 
the photon-magnon vertex diagrams in Fig. |^ is a posteriori verified. 

In order to take magnon-magnon interactions into account, we evaluate the scattering 

^1. 



vertex V'^ (see Fig. H) in the static limit for U ^ t. In Big geometry the only relevant 



contribution to is given by 

Vijqi,q2) = -l2J^,^^.,^, . (31) 

We decompose 7qi-g2 as 

791-92 = 7gi7g2 + lili + 7^iSf + 7^7^ , (32) 



where 



t = iisin(g^) + sin(gy)] and -ff = ^[sin(g^) - sin(gy)] . (33) 



7 

Due to the specific momentum dependence of the relevant photon-magnon vertex, Eq. (^8|), 
only the second term in Eq. ( ^21) contributes, because in performing the momentum sum 
in the vertex Eq. (^7]) all those parts of V^_^^ which are orthogonal to 7^ give a vanishing 
contribution (I]g7^7g = J2gl^l^^ = Hql'^jl^ = 0). Therefore, the effective remaining 
magnon-magnon interaction V^^j in Big geometry is given by 
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V,%{q,,q,) = -l2J^l^% . (34) 



Introducing the function 



1 1 

W^t) = II T7 II T^Xsc g", i^p)Xsci-Q, -q'^ - i^ph^' , (35) 

and solving Eqs. (pB]) and ( pT]) with V^"^ = Vfjj-{q[,q2) the Raman amphtude takes the 
form 

Finally, applying the cut prescription Eq. PT| ) yields the result for the non-resonant scat- 
tering intensity in Big geometry 

= | l-12JL(Ac. + .^) l • (3^) 

A comment is in order regarding the calculation of L{iujt). We note that the denominator 
of xtci^'^^) -^l- (HID contains two poles at u; = -^Ugyj^q) corresponding to forward and 
backward propagating spin wave excitations. The vertices in Fig. ^, however, account 
only for magnon number conserving scattering processes. For a consistent evaluation of the 
Raman amplitude it is required to retain only the forward propagating parts of the dynamic 
susceptibilities in Eq. ( pSf ) for L{iujt) p5| . 

With this restriction for the evaluation of L{iujt) the resulting Big scattering intensity is 
shown in Fig. |^ (solid line). 

As expected, the inclusion of magnon-magnon interactions removes the divergency of 
the intensity for zone boundary magnon pairs and leads to a single almost symmetric two- 
magnon peak around Au; = 2.6J. For comparison we have also shown the result of the 
Loudon-Fleury theory in Fig. |^ obtained from a spin wave theory analysis to order (1/5*)^. 
We note that an analogous but much more tedious expansion in the framework of the 
half-filled Hubbard model has been performed in Ref. pGf by extending the model to 25" 
equivalent orbitals on each site. 
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0.0 ' ■ ' — — ' ■ ' ■ 

0.0 1.0 2.0 3.0 4.0 

ACO [J] 

FIG. 5. Raman intensity in arbitrary units in the non-resonant case in the strong coupHng 
hmit with (sohd hne) and without (long dashed hne) final state magnon-magnon interactions. The 
short dashed line shows the result of the Loudon-Fleury theory with magnon-magnon interactions 
included. The transferred photon frequency Aa; is given in units oi J = 4:f/U. 

The line shape of the obtained spectra in Fig. ^is apparently distinctly different from the 
experimental data. As we argue in this paper the reason for this discrepancy is due to the 
neglect of resonance phenomena in the photon-magnon vertex function V; the non-resonant 
limit (2A — tUj) ^ J for which the Loudon-Fleury theory is valid is not the relevant situation 
for the experimental Raman data on cuprate antiferromagnets. 

VII. RESONANT RAMAN SCATTERING 

In this chapter we consider the resonant case when the incoming photon frequency Ui 
is comparable to the SDW energy gap to within a typical magnon energy, [ui — 2A) ~ J. 
In this regime strong resonant enhancement of the vertex function becomes important due 
to the combination of photon induced interband transitions and the creation of a magnon 
pair. Resonances appear in all four vertex diagrams in Fig. |^ which contribute to the total 
vertex function V. By inspection we find the strongest resonance to arise from a scattering 
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process that is contained in the rhomb-shaped vertex diagram A (Fig. ^). Alternatively to 
Eq. ( |Al| ) this contribution to the vertex function is more conveniently expressed in terms 
of SDW conduction and valence band quasi particle propagator matrices which allows more 
easily to identify the physical process responsible for the strongest resonance. Explicitly we 
rewrite the algebraic result for the diagram A as 

11'^ 
VA{q, iuip., iu!t, i^p) = J^Yl ^A^, g)Tr [ G^{k; iQ^) 

p(fc, k)G^{k, iVLn + i^iJi)v{k] k - q)G"'{k - q; iVtn + ioJ^ - iuOp) 

p{k - q, k - q)G'^ {k - q; iiln + i^^t - i^^p)pik - q, k) ] , (38) 

where we have introduced the coherence factor matrix 



piKk') 



V ^k,k' ~^k,k' ) 



(39) 



with matrix elements 



The symmetry factor M^ik^q) depends on the scattering geometry and is defined in Eq. 
( |A5| ) in the appendix. Carrying out the required matrix multiplications in Eq. (|38D leads to 
16 combinations of conduction {G'^) and valence band propagators {G""). Among them we 
single out the term which leads to the by far strongest resonant enhancement, i.e. the term 
containing the product G'"G'^G'^G'" , and neglect in the following all other combinations with 
weaker resonances. The physical process underlying the dominant resonance term is shown 
in Fig. 1^. The incoming photon creates a particle-hole pair by exciting an electron from the 
valence into the conduction band. This initial excitation then decays into a particle-hole 
pair with lower energy by the creation of two magnons with zero total momentum. Finally, 
the particle and the hole recombine under emission of the outgoing photon. 
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W\f\Pi q 



ca 



f\AAM-q 



FIG. 6. Scattering process which leads to the strongest resonant enhancement of the vertex 
function. The dashed Unes denote the incoming and outgoing photons, the wiggly hnes denote the 
magnons. E'"{k) and E^{k) mark the valence and conduction bands, respectively. 1.) Excitation of 
a particle- hole pair created by absorption of the incoming photon with frequency oJi. 2.) Creation 
of a magnon pair. 3.) Recombination of the particle-hole pair and emission of the outgoing photon 
with frequency ojf. 

After analytic continuation the vertex function contribution of the resonant scattering 
process shown in Fig. is given by 



m 



k,k 



{coi - 2E{k) + i5) 



X 



^k,k-q ^k-q,k ^k-q,k-q 



- E{k) -E{k-q)+ i5){uji - Alu - 2E{k - q) + i5) 



(41) 



Here we have already identified Au = uoi—uof = 2ujsw{<l) enforced by the energy conserving 
5-functions contained in the imaginary parts of the transverse dynamic susceptibilities (see 
Eq. (|T6|) ). The strong resonant enhancement of the vertex part results from the possible 
simultaneous vanishing of its three energy denominators and has therefore been termed a 



"triple resonance" in Ref. [|T5l. A triple resonance occurs if the following conditions hold: 



Ui = 2E{k) 



uji - Auj = 2E{k - q) . 



(42) 



The triple resonance conditions have been studied in detail analytically in Ref. |T5| by 
Chubukov and Frenkel. We have confirmed their results by solving numerically the triple 
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resonance equations for uji and Acj for a sequence of magnon momenta g in the MBZ. The 
numerical results are shown in Fig. |^ together with the position of the two-magnon peak 
in the non-resonant case. The boundary of the triple resonance region, i.e. the region in 
the (cui, Aa;) plane, where the resonance conditions can be matched, is determined by the 
high symmetry directions in the MBZ. Each point in the figure corresponds to a magnon 
momentum q and marks the frequencies Acj, uji for which the triple resonance equations 
(^) can be solved for electronic momenta k in the MBZ. The vertex function diverges 
for those values of Ui and Acu, which lie in the indicated region (Fig. ^ (as long as quasi 
particle lifetime effects are neglected). 
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FIG. 7. Numerical solution of the triple resonance conditions. Each point in the figure cor- 
responds to a pair {oji — 2A, Alij) where momenta k and q exist so that the resonance conditions 
are simultaneously fulfilled. The solid line marks the position of the two-magnon peak in the 
non-resonant case. 

In our slightly different calculational scheme we have so far reconfirmed the results of 
Chubukov and Frenkel for the location of the divergency of the vertex function arising from 
the vanishing of the denominators in Eq. (0). The obvious important task is now to explore 
the change of the two-magnon Raman intensity profile due to the triple resonance vertex 
function. 
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In the non-resonant case we essentially exploited the 7^ momentum dependence of the 
photon-magnon vertex function in the strong coupling limit for the calculation of the scat- 
tering intensity in Big geometry. In this case the 7^ momentum dependence arises simply 
from the Big symmetry factor MA{k,q) and from neglecting the momentum dependence 
of the SDW quasi particle propagators to leading order in t/U. In the present case the 
symmetry factor is left unchanged but the structure of the resonant vertex function mod- 
ifies the momentum dependence of VA{q)- Nevertheless, numerical calculations show that 
for intermediate photon frequencies 3 < (cuj — 2 A) / J < 5 the g^dependence of in Big 
geometry still follows very closely the 7^ form. In order to solve the Bethe-Salpeter equation 
for the renormalized vertex F even with a triple resonant photon-magnon vertex function 
we project on the 7|-channel, i.e. we single out the 7^ symmetry component of by 
introducing 

Vfiq, a;,, Au) = g{u„ Au) , (43) 

where 

E^^7|Vl^(g,^.,A^) .... 
g{uJi,Auj) = — -2 . (44) 

Using the projection of the triple resonance vertex term instead of the calculation 
of the Big scattering intensity proceeds completely analogous to chapter VI and leads to 

,45) 

The result for the resonant Big Raman intensity thus factorizes into the absolute square of 
the 7^ symmetry component of the vertex function and the two-magnon part which remains 
unchanged. The spectrum therefore consists of two separate contributions of distinct origin: 
the two-magnon peak at Auj ^ 2.6J and the triple resonance peak which appears well 
above the two-magnon peak close below 4J for an intermediate photon frequency range 
3 < (cui - 2A)/J < 5 (see Fig. 0). 

Clearly, the magnitude of the triple resonance structure in the two-magnon intensity 
profile depends on the strength of the quasi particle damping resulting from self-energy 
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corrections due to residual interactions between the SDW quasi particles ||2^. Similarly, 



its precise location depends on the renormalized band dispersion. Here, for the purpose of 
demonstrating the consequence of the triple resonance vertex on the two-magnon line shape, 
we model the effect of quasi particle damping by adding a finite imaginary part to the energy 
denominators of the triple resonant vertex function. For the choice of a typical quasi particle 
lifetime we are guided by the results of a self-consistent non-crossing calculation of the self- 
energy correction in the SDW state of the half-filled Hubbard model ||2^ . The result for the 
two-magnon intensity in Big geometry with magnon-magnon interactions and a broadened 
triple resonance vertex function is shown in Fig. 

VIII. DISCUSSION 

Fig. ^ shows a comparison of our result for the two-magnon Raman intensity in Big 
geometry to the experimental spectrum for La2CuOii taken from Ref. In order to isolate 
the two-magnon signal a background was substracted from the experimental data as in 



Ref. ||3T|. The calculated lineshape is dominated by the two-magnon peak but its high 
energy shoulder results solely from the triple resonance enhancement in the photon-magnon 
vertex function. The charge transfer energy gap in La2CuOi is about 2eV as deduced from 
measurements of the optical conductivity p8[. So for the laser photon frequency Ui = 2.55eV 



used in the Raman experiment we have cjj— 2A ~ 550meV. For the calculation of the Raman 
intensity we chose J = 1200cm~^ such that the two-magnon peak frequency coincides with 
the experimental value and tUj — 2A = 3. 6 J ~ 550meV. For intermediate photon frequencies 
{uji — 2A) ~ 4J we obtain a lineshape that is in fair agreement with the experimental 
spectra not only for La2Cu04 but also for the AF double-layer compound YBa2Cu30Q,i 
near resonance as shown in Fig. p|b for a series of photon frequencies [E^ . 
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FIG. 8. a) Raman scattering intensity in arbitrary units in Big geometry as calculated from 
Eq. (^) with (uji — 2A) = 3.6J for a lattice with 100 x 100 sites. For the two-magnon term an 
imaginary broadening of i5 = iO.OdJ was used in the frequency denominators, while the vertex part 
was evaluated with 5 = 0.4J. The jagged line displays the experimental spectrum of La2Cu04, 



taken at room temperature with a laser frequency of = 2.55eV Q. As in Ref. |31| a background 
was substracted from the Raman data. For comparison with the data the magnetic energy scale 
was set to J = 1200cm~^ and the peak intensities were scaled to coincide, b) Experimental spectra 
of YBa2Cu30Q,i at T = 5K for various laser energies as indicated in the figure (reproduced from 
Ref. M). 
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Besides the high energy shoulder other experimental features in the scattering intensity 
find a natural explanation as well by the interplay between the two-magnon peak and the 
triple resonance: 

(a) Two separate resonance frequencies. Recent measurements of the Raman intensity by 
Blumberg et al. |^ show an enhancement of the absolute Raman intensity for two different, 
well separated photon frequencies Ui. In our analysis the Raman intensity as a function of 
Ui is expected to increase if the frequencies for the two-magnon and the triple resonance 
peak merge. This does, in fact, happen for two distinct frequencies {ui — 2A) ^ 2.6J and 
[ui — 2 A) ^ 8 J. These frequencies can be read off from Fig. |^ from the intersections of the 
two-magnon peak frequency (horizontal line) with the triple resonance region. 
(6) uj'^'^^ > 2A. The same arguments as in (a) apply also for understanding the experimental 
observation that resonance occurs for incoming photon frequencies well above the insulating 



energy gap 2 A. This has already been pointed out previously in Ref. [ p9[] . 
(c) Changing line shape. The experimental spectra in Fig. show a decreasing asymme- 
try with increasing photon frequency tUj. The two spectra for YBa2Cu30Q,i taken at the 
lowest frequencies contain two well separated contributions, whereas the total line shape 
of the spectrum of the largest laser frequency is nearly symmetric and the two peaks are 
indistinguishable. (For a detailed discussion of these spectra see Ref. |^.) Within our 



calculation, the almost symmetric intensity profile corresponds to the near coincidence of 
the triple resonance and the two-magnon peak frequencies. 

Fig. H a shows Raman data of La2CuOi, which were taken at room temperature [Q. The 
most noticable difference to the calculated intensity is the larger peak width in the experi- 
mental spectrum. There are several possible explanations for this discrepancy. First of all 
one might wonder about the influence of flnite temperatures, since our calculation was per- 
formed at T = 0. However, a comparison between the low and high temperature (T ~ T/vee«) 
spectra in Fig. ||a and Fig. ^b for La2CuOi and Y Ba2Cu?jOQ,i, respectively, reveals that 
the thermal effects are of minor importance since the spectra have a comparable peak width. 
Clearly, lifetime effects due to magnon-phonon and also magnon-magnon interactions in par- 
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ticular for the zone boundary magnons will broaden the two-magnon Raman signal. The 
Raman spectrum in Fig. ^ja was evaluated on a 100 x 100 lattice and a small broadening 
of i5 = i0.09J in the energy denominator of the dynamic spin susceptibility was introduced 
only for calculational purposes. The effects of magnon-phonon and magnon-magnon inter- 
actions were neglected in the present work and the sharper peak of the calculated Raman 
intensity in comparison to the data is expected on physical grounds. The intention of our 
present analysis was, instead, to show that the resonant frequency dependence of the vertex 
function for the photon-magnon coupling - similar to and in agreement with the recent work 
by Chubukov and Frenkel JTSf - leads to a qualitative change of the two-magnon line shape 



and provides a natural explanation for the resonance phenomena observed in light scattering 
experiments on cuprate antiferromagnets. 

Furthermore we neglected the contribution from four-magnon and higher order scattering 
processes. Canali and Girvin have attempted to include the four-magnon contribution in 



the framework of the Heisenberg model and the Loudon-Fleury coupling Hamiltonian pi 



Although the first three moments of the resulting scattering intensity are in good agreement 



with experiment and also with the series expansion results by Singh et al. PT[, the exper- 
imental line shape is not successfully explained. The integrated four-magnon intensity is 
only 2.9% of the two-magnon intensity, and if this estimate is correct, the four-magnon con- 
tribution will barely affect the overall Raman intensity profile. Furthermore, the Canali and 
Girvin four-magnon peak frequency is 2.5 times larger than the two-magnon peak frequency 
and therefore one has to conclude that four-magnon scattering cannot be responsible for the 
structure in the Raman intensity near 4000cm~^. 

The effects of phonon-magnon interaction have been studied previously in Refs. [p^-|3^ 
for Raman scattering in the Heisenberg antiferromagnet. As expected, the damping of 
zone boundary magnons does indeed lead to a significant broadening of the two-magnon 
signal. A qualitatively smaller contribution to the damping arises also from magnon-magnon 



interactions Nori et al. considered the effect of a random Gaussian variation SJij of 



the exchange couphng Jq assumed to originate from a distortion of the crystal lattice due to 
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low frequency lattice vibrations ||3J]. Using numerical techniques Nori et al. showed that in 
this model a broad and asymmetric two-magnon line shape with an enhancement of spectral 
weight at higher energies emerges if the mean deviation < 5Jij > is as large as Jo/2. Such 
an enormous variation of the exchange coupling, however, appears physically unreasonable. 
It is far from being clear whether more moderate values of < 6Jij > are sufficient for this 
mechanism to explain the line shape of the two-magnon Raman spectra. We emphasize that 
in all these works p2H3^ the resonant nature of the two-magnon scattering process was not 
taken into account. 



IX. SUMMARY AND CONCLUSIONS 

Based on the single-band Hubbard model at half filling we have performed a microscopic 
analysis of two-magnon Raman scattering in a SDW antiferromagnet. In a diagrammatic 
formulation we have explicitly taken into account the structure of the photon-magnon cou- 
pling vertex. This allows to explore both, the non-resonant [ui ^ 2A) and the resonant 
case {uji — 2 A ~ J) for the light scattering intensity, where the latter is the relevant limit for 
Raman experiments on undoped cuprate antiferromagnets. For the non-resonant case in the 
strong coupling limit the results of the conventional Loudon-Fleury theory for two-magnon 
Raman scattering in the Heisenberg antiferromagnet are almost quantitatively reproduced 
verifying our selection of the relevant vertex diagrams. In the resonant regime we identified 
the physical process which yields the strongest diverging part of the photon-magnon vertex 
function and analyzed the conditions and consequences of a triple resonance. For photon 
frequencies well above the SDW energy gap triple resonance occurs for transferred photon 
frequencies larger than the two-magnon peak frequency leading to a high energy shoulder of 
the two-magnon intensity profile. These results confirm the conclusions of the triple reso- 
nance theory by Chubukov and Frenkel for resonant two-magnon Raman scattering. In 
this theory effective electron-magnon vertices were constructed from the requirement that 
these vertices lead to the strong coupling RPA form of the transverse dynamic spin sus- 
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ceptibility. In contrast to the work of Chubukov and Frenkel we have performed a more 
systematic microscopic derivation for the two-magnon scattering intensity and we have eval- 
uated the intensity profile by explicitly taking into account the resonant structure of the 
photon-magnon coupling vertex. 

By using unrenormalized SDW quasi particle propagators the triple resonance leads 
to a true divergence of the vertex function. Self-energy corrections, however, will remove 
the divergence leaving a finite enhancement in the Raman intensity at the triple resonance 
frequency. These self-energy corrections have been calculated in previous work We have 
here for simplicity modelled quasi particle lifetimes effects by adding a constant damping 
term to the propagators in the triple resonance vertex function. A more reliable quantitative 
estimate for the strength of the triple resonant enhancement will require a refined treatment 
with the inclusion of self-energy corrections. 

In conclusion, we found that the resonant frequency dependence of the photon-magnon 
vertex function gives rise to an enhancement of the high energy side of the two-magnon 
Raman peak in Big scattering geometry. The combination of resonant transitions between 
the SDW quasi particle bands and magnon pair excitations provides a microscopic basis for 
understanding the resonant Raman scattering experiments on cuprate antiferromagnets. 
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APPENDIX: 



For the calculation of the vertex diagrams shown in Fig. |^ we use the bare SDW sin- 
gle particle propagator Gg (in c-fermion representation) in 2 x 2 matrix notation. In the 
finite temperature Matsubara formalism the vertex functions of the four diagrams are given 
explicitly by 



I E ^E'^a(^, q)Tr [G^(fc; tQ,,)(T, G^(fc; zfi„ + too,) 



if!,i,jf!„ kl 



cr^Gg '"(A; - q; iQn + ii^t] 



UTi 



G^(/i«a 



Gq (Z; ifln + 2Co't)GQ °"(/ — q; ifln + iuJt — ioJp 



(A2) 



Vc{q, iuJ^, iut, iUp) = -5 E T?E ^i?(^)Tr Gg (A;; iil„)cr^ Gg (/c; iQ^ + iuJ^)fT, 



N 

Gq(A;; iQn + ^i^i)Go'^(fc — q; iQn + ^^^t — i^p 



(A3) 



^ E ^E'^i?(^)Tr [G^(^; zfi„)o-, G^(^; tQ^ - iuj,)(t, 
Gq (fc; 2fi„ + i(x't)Go '^(/c — (f; + icUf — iup 



(A4) 



The symmetry factors of the basic scattering vertices from the coupling of the photon vector 
potential to the electron current density have been combined into the functions 



MA{k,q} 
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.yUJiujJ dka dkp 



(A5) 



MB{k) = MA{k,Q) 



(A6) 



Note that the factor de^jdka changes sign when umklapp scattering (A; — /c + Q) has taken 
place along the fermion lines. In the 2x2 matrix formulation the alteration of the sign is 
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conveniently taken into account by inserting the Pauli matrix (T^. Including all diagrams 
with reversed direction of the fermion lines yields a factor two in the total vertex function. 
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